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Recent theoretical progress in understanding high-energy scattering beyond the mean 
field approximation is reviewed. The role of Lorentz invariance and pomeron loops in 
the evolution, the relation between high-energy QCD and statistical physics and results 
for the saturation momentum and the scattering amplitude are discussed. 

1 Introduction 

The high-energy scattering of a dipole off a nucleus/hadron in the mean field approximation 
is decribed by the BK-equation pQ . The main results following from the BK-equation are the 
geometric scaling behaviour of the scattering amplitude and the roughly powerlike energy 
dependence of the saturation scale which are both nicely supported by the HERA data. 

The recent progress consists in understanding small- a; dynamics (near the unitarity limit) 
beyond the mean field approximation, i.e., beyond the BK-equation. A first step beyond the 
mean field approximation was done in [5] where the BFKL evolution in the scattering process 
was enforced to satisfy natural requirements as unitarity limits and Lorentz invariance. The 
result was a correction to the saturation scale and the breaking of the geometric scaling at 
high energies. Afterwards a relation between high-energy QCD and statistical physics was 
found [3] which has clarified the physical picture of, and the way to deal with, the dynamics 
beyond the BK-equation. It has been understood that gluon number fluctuations from one 
scattering event to another and the discreteness of gluon numbers, both ignored in the 
BK evolution and also in the Balitsky-JIMWLK equations [I], lead to the breaking of the 
geometric scaling and to the correction to the saturation scale, respectivelly. New evolution 
equations [6j [7] , which describe Pomeron loops, have been proposed to account for the 
above effects. Very recently possible effects of Pomeron loops on various observables piUlfTU] 
have been studied in case they become important in the range of collider energies. 

In the following sections I will show the recent developments in some detail by considering 
equations and results in and beyond the mean field approximation. 

2 Mean field approximation 

Consider the high-energy scattering of a dipole of transverse size r off a target (hadron, 
nucleus) at rapidity Y = ln(l/x). The rapidity dependence of the T-matrix in the mean field 
approximation is given by the BK-equation which has the schematic structure (transverse 
dimensions are suppressed) 

d Y T = a s [T-TT] . (1) 

The linear part of the BK equation, i.e., the BFKL equation, gives the growth of T with 
rapidity Y whereas the non-linear term TT tames the growth of T in such a way that the 
unitarity limit T < 1 is satisfied. 
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One of the main results following from the BK-equation is the geometric scaling behaviour 
of the T-matrix [UJ 

T(r 7 Y)=T(r 2 Q 2 s (Y)) , (2) 

where Q S (Y) is the so-called saturation momentum defined such that T(r ~ 1/Q S ,Y) be of 
0(1). Eq. ([2]) means that the T-matrix scales with a single quantity r 2 Q 2 (Y) rather than 
depending on r and Y separatelly. This behaviour implies a similar scaling for the DIS cross 
section, a^* p {Y, Q 2 ) = <t^f(Q 2 /Q 2 s {Y)), which is supported by the HERA data. 

Another important result that can be extracted from the BK-equation is the rapidity 
dependence of the saturation momentum (leading- Y contribution) |12| , 



Q 2 (Y) = Q 2 Exp 



2a s N c x(Ap) y 

7T 1 - An 



(3) 



where x(A) is the BFKL kernel and A = 0.372. 

The shape of the T-matrix resulting from the BK-equation is preserved in the transition 
region from weak to strong scattering, < T < 1, with rising Y (front of the travelling 
wave): The saturation region at r > 1/Q S (Y) where T ~ 1 however widens up, including 
smaller and smaller dipoles, due to the growth of the saturation momentum. The situation 
changes, as we will, once gluon number fluctuations are taken into account. 



3 Beyond the mean field approximation 

3.0.1 Lorentz invariance and unitarity requirements 

Let's start with an elementary dipole of size r\ at rapidity y = and evolve it using the 
BFKL evolution up to y = Y. The number density of dipoles of size r 2 at Y in this dipole, 
n { r i, r 2,Y), obeys a completeness relation 

n( ri ,r 2 ,Y) = [ n(n,r,Y/2) n(r,r 2 ,Y/2) (4) 

where on the right hand side the rapidity evolution is separated in two successive steps, 
y = ^ y = Y/2 ^ y = Y. With 

T(n,r 2 ,Y) ~ c a 2 r\ n(r 1 ,r 2 ,Y) (5) 

eq.© can be approximately rewritten in terms of the T-matrix as 

T( ri ,r 2 ,Y)) *^fdp (1 r(n,r,y/2)) (L T(r,, 2 ,y/2)) (6) 

where p = ln(rQ/r 2 ). In Ref. [5] it was realized that the above completeness relations, 
or, equivalently, the Lorentz invariance, is satisfied by the BK evolution only by violating 
unitarity limits. This can be illustrated as follows: Suppose that r 2 is close to the saturation 
line, r 2 ~ l/Q s (Y/2), so that the left hand side of Eq.© is large. On the right hand side 
of Eq.© it turns out that T(ri,r,Y/2)/r 2 is typically very small in the region of p which 
dominates the integral. This means that T(r, Y/2)/r% must be typically very large and 
must violate unitarity, T(r, r 2 , Y/2) ^> 1, in order to satisfy ©. 
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The simple procedure used in Ref. [5] to solve the above problem was to limit the region 
of the p-integration in Eq.® by a boundary p2(Y/2) so that T(r, r%, F/2)/r| would never 
violate unitarity, or T(ji , r, Y/2) / r 2 would always be larger than a 2 . The main consequence 
of this procedure, i.e., BK evolution plus boundary correcting it in the weak scattering 
region, is the following scaling behaviour of the T-matrix near the unitarity limit 

T(r Y) = T ( l ^Q 2 s(Y)) ) (7) 
and the following energy dependence of the saturation momentum 



Q 2 S (Y) = Ql Exp 

with 



2a s iV c x(A ) v L ^x"(Ao) 



-Y 1 



1 - A V 2(Ap) 2 x (A ] 



(8) 



113 1 

Ap = — In — + — In In — + const. . (9) 

1 - A aj 1 - A aj. 

Eq.® shows a breaking of the geometric scaling which was the hallmark of the BK equa- 
tion (cf. Eq.([2|)) and Eq.© shows the correction to the saturation momentum due to the 
evolution beyond the mean field approximation (cf. Eq.((3])). 



3.0.2 Statistical physics - high density QCD correspondence 

The high energy evolution can be viewed also in another way which is inspired by dynamics of 
reaction-diffusion processes in statistical physics [3J. To show it, let's consider an elementary 
target dipole of size r\ which evolves from y = up to y = Y and is then probed by an 
elementary dipole of size r, giving the amplitude T(ri,r, Y). It has become clear that the 
evolution of the target dipole is stochastic leading to random dipole number realizations 
inside the target dipole at Y, corresponding to different events in an experiment. The 
physical amplitude, T(ri,r, Y), is then given by averaging over all possible dipole number 
realizations/events, T(ri,r, Y) = (T(ri,r, Y)), where T(r\,r, Y) is the amplitude for dipole 
r scattering off a particular realization of the evolved target dipole at Y . 

The mean field descripton breaks down at low target dipole occupancy due to the discrete- 
ness and the fluctuations of dipole numbers. Because of discreteness the dipole occupancy 
can not be less than one for any dipole size. Taking this fact into account by using the 
BK equation with a cutoff when T becomes of order a 2 s [3J , or the occupancy of order one 
(see Eq.®), leads exactly to the same correction for the saturation momentum as given in 
Eq.(j3j. The latter cutoff is essentially the same as, and gives a natural explanation of, the 
boundary used in Ref. [2) and briefly explained in the previous section. 

The dipole number fluctuations in the low dipole occupancy region result in fluctuations 
of the saturation momentum from event to event, with the strength 

a 2 = (pi) - ( Ps ) 2 = const. (10) 

extracted from numerical simulations of statistical models. The averaging over all events 
with random saturation momenta, in order to get the physical amplitude, causes the breaking 
of the geometric scaling and replaces it by the scaling law 

(T(r,Y))=T(^mi) . (11) 
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This equation differs from Eq.0 since Eq.0 misses dipole number fluctuations. 



3.0.3 Pomeron loop equations 

It was always clear that the BK equation does not include fluctuations. However, it took 
some time to realize that also the Balitsky-JIMWLK equations do miss them. As soon as 
this became clear (first Ref . in [5] ) , the so-called Pomeron loop equations |Hj have been 
constructed, aiming at a description of fluctuations. They can be written (schematic way, 
transverse dimensions ignored) as a stochastic equation of Langcvin-type, 



d Y T = 



T - TT + a, VT v 



(12) 



or as a hierarchy of coupled equations of averaged amplitudes, the first two of them reading 

d Y (T) = a s [(T)-(TT}} 

dy (TT) = a s [(TT) - (TTT) + a 2 s (T)] . (13) 

The last term in Eq. p2[) . containing a non-Gaussian noise is, is new as compared with the 
BK equation and accounts for the fluctuations in the dipole number. Eq. (|13p reduces to the 
BK equation only in the mean field approximation, i.e., if (TT) = (T)(T). The hierarchy in 
Eq. (fT3")) . as compared with the Balitsky-JIMWLK hierarchy, involves in addition to linear 
BFKL evolution and pomeron mcrgings, also pomeron splittings, and therefore pomeron 
loops. The three pieces of evolution are represented by the three terms in the second equation 
in Eq. (fT3|) . respcctivclly, in the case where two dipoles scatter off a target. 

It isn't yet clear at which energy fluctuation effects start becoming important. The 
results shown in the previous sections, Eq.® and Eq. (|lip , are valid at asymptotic energies. 
A solution to the evolution equations, which is not yet available because of their complexity, 
would have helped to better understand the subasymptotics. However, using the methods 
outlined in the previous subsections, phcnomcnological consequences of fluctuations in the 
fixed coupling case have been studied, for example for DIS and diffractivc cross sections [5], 
forward gluon production in hadron-hadron collisions [lOj and for the nuclear modification 
factor R p A [5], in case fluctuations become important in the range of LHC energies. 
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